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(57) ABSTRACT 

A system with applications in pattern recognition, or classi- 
fication, of DNA assay samples. Because DNA reference and 
sample material in wells of an assay may be caused to fluo- 
resce depending upon dye added to the material, the resulting 
light may be imaged onto an embodiment comprising an 
array of photodetectors and an adaptive neural network, with 
applications to DNA analysis. Other embodiments are 
described and claimed. 
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BIOCHIP MICROSYSTEM FOR 
BIOINFORMATICS RECOGNITION AND 
ANALYSIS 

PRIORITY CLAIM 5 

This application claims the benefit of U.S. Provisional 
Application No. 60/856,512, filed 3 Nov. 2007. 

GOVERNMENT INTEREST 10 

The invention described herein was made in the perfor- 
mance of work under a NASA contract, and is subject to the 
provisions of Public Law 96-517 (35 USC 202) in which the 
Contractor has elected to retain title. 

FIELD 

Embodiments relate to recognition of DNA material using 2Q 
an optical and electronic system. 

BACKGROUND 

Genetic material may be analyzed by placing DNA (Deox- 2 5 
yribonucleic acid) material in an array of wells (dots). Poly- 
merase chain reaction (PCR) amplification is often used in 
genetic analysis, where the PCR amplification augments the 
amount of DNA material placed in a well. Fluorescent dyes, 
such as CY3 or CY5, may be added to the DNA material, so 30 
that it fluoresces when excited by monochromatic light. 
Because with PCR amplification there may be different 
growth rates of DNA material from well to well, sample and 
reference channels may be set up whereby in each well, there 
is reference DNA and sample DNA. Fluorescent dye of one 35 
type may be used for the sample DNA, and fluorescent dye of 
another type may be used for the reference DNA. 

This method also reduces the sources of variability and 
noise due to various aspects of an individual spot that affect 
both specimens (DNA sample and reference) similarly. In 40 
order to accurately calculate the density of the sample DNA 
material in a particular well after PCR amplification, the 
integral of the total fluorescence intensity (presumably rep- 
resenting the density of the DNA material inside the well) 45 
from the topological profile of the well is usually computed. 
The logarithmic value of the ratio of the two intensities of the 
fluorescent dye labeled specimens (one value for the sample 
specimen, the other value for the reference specimen) mea- 
sured from the same well is calculated based on the assay’s 50 
fluorescence image. The ratio of the two intensities would 
provide the normalized population of the gene material in the 
well, disregarding the initial population density. 

In most of the available commercial solutions, the assay’s 
fluorescence image is usually scanned by a color scanner with 55 
high resolution and then transferred to a computer for image 
analysis. The profile analysis software usually computes the 
normalized intensity of each well sequentially. The intensity 
of the fluorescence is usually relatively low. Using higher 
excitation light intensity or increasing detection time may 60 
lead to brighter fluorescence patterns. However, lower power 
consumption and faster detection may be preferable. Further- 
more, some fixed-pattern noises in the input pattern may exist 
(e.g., fixed pattern noises created by scattered lights, or non- 
uniformity of the detector array response). These noises may 65 
introduce errors in the measurement of the density of the 
DNA materials 
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The development of low-cost portable instruments for rap- 
idly analyzing genetic assays in noisy environments and with 
relatively low intensity of fluorescence would be of utility in 
medical services. 

BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 illustrates a system for recognizing a pattern among 
a DNA sample assay according to an embodiment. 

FIGS. 2 and 3 illustrate adaptive neural networks accord- 
ing to an embodiment. 

FIG. 4 illustrates functional blocks of a VLSI implemen- 
tation according to an embodiment. 

FIG. 5 illustrates a functional block diagram of a pixel 
according to an embodiment. 

DESCRIPTION OF EMBODIMENTS 

In the description that follows, the scope of the term “some 
embodiments” is not to be so limited as to mean more than 
one embodiment, but rather, the scope may include one 
embodiment, more than one embodiment, or perhaps all 
embodiments. 

FIG. 1 illustrates various components of an embodiment in 
an exploded view. Assay array 102 comprises an array of 
dual-labeled gene wells (dots). For simplicity of illustration, 
only one well, labeled 104, is illustrated. Each well comprises 
reference DNA material and sample DNA material. The 
method of PCR amplification is applied to the array of wells 
in which a first dye is added to the sample DNA material and 
a second dye is added to the reference DNA material. When 
excited by light, the reference and sample DNA material may 
be made to fluoresce. The intensity of the fluorescence is 
indicative of the amount of DNA sample or reference mate- 
rial. The excitation light for the sample DNA material will for 
many instances have a different spectrum than the excitation 
light for the reference DNA material. The light given off by 
the sample DNA material will usually have a different spec- 
trum than the light given off by the reference DNA material, 
depending upon the dyes, and usually the light given off by 
the reference and sample DNA material will have a spectrum 
different from their respective excitations. For some embodi- 
ments, the excitation may be monochromatic light. The 
monochromatic light may be realized by using optical notch 
filters in front of a relatively broad light source. 

For example, for some embodiments CY3 and CY5 dyes 
may be used, where the DNA material to which CY3 has been 
added is excited by monochromatic light having a center 
frequency (wavelength) of 535 nm and a bandwidth of 1 0 nm, 
and the DNA material to which CY 5 has been added is excited 
by monochromatic light at 625 nm with a bandwidth of 10 
nm. For such embodiments using CY3 and CY5 dyes, when 
fluorescing the CY3 dye gives off green light having a peak 
value at 570 nm, and the CY5 dye gives off red light having a 
peak value at 670 nm. The excitation of wells with CY3 and 
CY 5 may be performed concurrently. Assay array 1 02 may be 
front-side illuminated, or backside illuminated if transparent, 
for example. 

The term monochromatic is a term of art, where of course 
in theory no excitation source is purely monochromatic. In 
the examples given, the bandwidth is less than about Vso of the 
center frequency. For some embodiments, the excitation need 
not be monochromatic in this sense. 

Lens system 106 images the light from the array of wells 
onto sensor array 110. (Lens system 106 may comprise more 
than one lens element.) Sensor array 110 comprises an array 
of pixels, but for simplicity of illustration only one pixel, 
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labeled 112, is shown in FIG. 1. Pixel 112 comprises two 
optical filters, 112a and 1126, where optical filter 112a has a 
pass band to allow the fluorescence from the DNA sample 
material to pass through but to substantially reject light out- 
side the frequency range of this fluorescence. Similarly, opti- 
cal filter 1126 has a pass band to allow the fluorescence from 
the DNA reference material to pass through but to substan- 
tially reject light outside the frequency range of this fluores- 
cence. 

For example, for embodiments using CY3 and CY 5 dyes as 
discussed above, one of the optical filters, say 112a, is a 
thin-film micro -optical-filter having a passband centered at 
about 580 nm with a bandwidth of 40 mil, and the other 
optical filter, say 1126, is a thin-film micro -optical -filter hav- 
ing a passband centered at about 765 nm with a bandwidth of 
40 nm. In this way, a sensor below optical filter 112a is 
responsive to the DNA material having the CY3 dye, and a 
sensor below optical filter 1126 is responsive to the DNA 
material having the CY5 dye. 

An exploded view of the sensors and circuit for pixel 112 of 
the embodiment of FIG. 1 is provided as circuit 114, which 
may be referred to as a differential logarithm circuit. In the 
embodiment of FIG. 1, photodetector 116a is below optical 
filter 112a, and photodetector 1166 is below optical filter 
1126. In the particular embodiment of FIG. 1, photodetector 
116a is responsive to the light imaged from the sample DNA 
material, and photodetector 1166 is responsive to the light 
imaged from the reference DNA material. In the particular 
embodiment of FIG. 1, photodetectors 116a and 1166 are 
NPN photodetectors, but other embodiments may use other 
types of photodetectors. 

Photodetector 116a is connected in series with two tran- 
sistors, nMOSFETs (n- Metal-Oxide-Field-Effect-Transis- 
tor) 118a and 120a. The drain-source currents for nMOS- 
FETs 118a and 120a are substantially equal to the current 
sourced by photodetector 116a. The current sourced by pho- 
todetector 116a is proportional to the amplitude of the inci- 
dent light. Transistors 118a and 120a are each diode-con- 
nected. When operating in their sub-threshold regions, their 
gate-to -source voltages are substantially proportional to the 
logarithm of the current sourced by photodetector 116a, 
which in turn is proportional to the logarithm of the amplitude 
of the light incident on photodetector 116a. Similar remarks 
apply to photodetector 1166, and transistors 1186 and 1206, 
but where photodetector 1166 is responsive to incident light 
from the reference material. 

Differential transistor pair 112a and 1126, resistors 124a 
and 1246, and tail current transistor 126 form a differential 
amplifier, where the input signals are the gate voltages of 
transistors 118a and 1186, and the output voltage is taken at 
the drain of one of the transistors in the differential transistor 
pair. For the particular embodiment of FIG. 1, the drain on 
transistor 1226 is taken as an output port, labeled as V OLrr 
Denoting the amplitude of the sample incident light by A s and 
the amplitude of the reference incident light by the A^, the 
output voltage may be written as 



where K is some proportionality factor. It should be noted that 
the above equation for the output voltage is only approximate, 
and does not serve as an exact expression of the output voltage 
in terms of A s and A^. 
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Other embodiments may employ circuits different from the 
particular circuit illustrated in FIG. 1. For example, a circuit 
complementary to circuit 114 may be realized, where pMOS- 
FETs are used instead of nMOSFETs in circuit 114, and PNP 
5 photodetectors are used instead of photodetectors 116a and 
1166. Furthermore, some embodiments may employ other 
types of transistors, such as bipolar transistors. As another 
example, whereas the embodiment for circuit 114 illustrates 
two diode-connected transistors in series with each photode- 
10 tector, other embodiments may use a different number of 
diode-connected transistors. Other embodiments may use 
other types of differential amplifiers in place of the differen- 
tial amplifier represented by transistor pair 112a and 1126, 
resistors 124a and 1246, and tail current transistor 126. For 
15 example, some embodiments may employ active devices in 
place of resistors 124a and 1246, where such active devices 
have a relatively wide range linear impedance response, or 
other embodiments may employ a different configuration of 
transistors to provide the tail current provided by transistor 
20 126. 

Accordingly, pixel 112 may be represented by the func- 
tional blocks indicated in FIG. 5. Referring to FIG. 5, optical 
filter 112a has a first passband to pass to optical detector 502a 
light that has been filtered by the first passband. The combi - 
25 nation of optical detector 502a and logarithmic circuit 504a 
provides a voltage to differential amplifier indicative of the 
logarithm of the light intensity provided to optical detector 
502a. Similar remarks apply to optical filter 1126, optical 
detector 5026, and logarithmic circuit 5046, but where optical 
30 filter 1126 has a second passband different tuned to a different 
frequency spectrum from that of the first passband. The out- 
put of differential amplifier 506 is indicative of the logarithm 
of the ratio of the light intensities provided to optical detectors 
502a and 5026. 

35 Given an array of voltage signals, each voltage signal 
indicative of the logarithm of the ratio of the sample light 
amplitude to the reference light amplitude for a particular pair 
of wells from the assay array, embodiments may use an adap- 
tive neural network to classify the voltage signals. Classifi- 
40 cation may be viewed as pattern recognition. Some embodi- 
ments may employ an adaptive neural network structure such 
as that illustrated in FIG. 2, where the layer of neurons 202 is 
a layer of input neurons and the layer of neurons 204 is a layer 
of output neurons. Input neurons 202 pass on their input to the 
45 next layer of neurons, neurons 204, where neurons 204 per- 
form processing on their input. 

Shown in FIG. 2 is an enlarged view of a neuron 204, 
indicating a summation function 2 using weights {w z , i=l, 
2, . . . , M}, and a transfer function A. The integer M denotes 
50 the number of inputs to neuron 204, which for the embodi- 
ment of FIG. 2 is the number of input neurons 202. Denoting 
the inputs to neuron 204 as {x., i=l, 2, . . . , M}, the summation 
function provides an intermediate term h where 

55 h=2 i=1 M w iXi . 

The weights depend upon the particular output neuron per- 
forming the summation, but to avoid multiple subscripts this 
dependency is not explicitly indicated. For some embodi- 
6Q ments, a bias term b may be used, where 

h=J. i=l M wpc r b, 

where the bias term depends upon the particular output neu- 
ron. 

65 The intermediate term h is fed as input to the transfer 
function A to provide an output y, where in general y=A(h). 
For some embodiments, A(h) may take the following form: 
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(1 + exp {-h)) 1 h < 2 
-a\n(B(6 - h)) -2<h<0 

A(h) = 

adn(/?((5 + h)) 0 < h < 2 

(1 + exp(-/i)) _1 2 <h 

where a, (3, and 6 are constants. As a particular example for 
some embodiments, the values a=0.050095635, p=1000, and 
6=0.01 have been used in experiments. Sometimes, the trans- 
fer function is referred to as an activation function, which is 
the motivation for using the notation A. 

The particular transfer function described above may be 
termed a sigmoid-logarithmic transfer function, or piece- 
wise sigmoid-logarithmic transfer function. This is to be dis- 
tinguished from the relatively common sigmoid transfer func- 
tion where A(h)=(l+exp(-h)) _1 for all h. 

For some embodiments employing digital processing to 
perform the adaptive neural network function, the input volt- 
age signals to input neurons 202 are quantized by one or more 
analog-to-digital converters so that input neurons 202 and 
output neurons 204 operate in the digital domain. However, 
other embodiments may be mixed-signal systems, where 
some processing functions are performed in the analog 
domain, and some processing functions are performed in the 
digital domain. Some embodiments may be realized in which 
almost all functions, or all functions, are performed in the 
analog domain. For example, the weighted summation per- 
formed by a neuron, as well as the transfer function, may be 
performed in the analog domain using analog multiplier cir- 
cuits and analog summation circuits. 

For some embodiments, the adaptive neural network may 
comprise more than one processing layer, so that there are 
hidden layers. For example, illustrated in FIG. 3 is an adaptive 
neural network with two input neurons (layer 301), a hidden 
layer with three neurons (layer 3 02), and an output layer with 
two neurons (layer 303), where neuron 302' is a neuron in 
hidden layer 302 with transfer function A, and neuron 303* is 
an output neuron in layer 303 with transfer function A. The 
adaptive neural network illustrated in FIG. 3 is a feedforward 
network because there are no feedback paths from a neuron in 
one layer to another neuron in a preceding layer. Some 
embodiments, for example, may include feedback paths from 
layer 303 to layer 302 to implement a recursive adaptive 
neural network. 

In FIG. 3, the inputs are represented by x 1 and x 2 . These are 
provided as inputs to the hidden layer 302, such as for 
example neuron 302'. For neuron 302', it is understood that 
the summation operation operates on two inputs (i.e., the 
inputs x 1 and x 2 ). The output of any one neuron in hidden 
layer 302 is provided as an input to all neurons in output layer 
303. The output of neuron 302' is denoted as x/, where the 
index i=l , 2, 3 corresponds to neurons 302 a, 3026, and 302c, 
respectively. This output is provided as one of three inputs to 
each neuron in output layer 303, such as neuron 303’. It is 
understood that the summation sign in neuron 303’ operates 
on three inputs (i.e., x 1 ', x 2 ', and x 3 '). For some embodiments, 
the transfer function A may be that as described in the previ- 
ously displayed expressions for A. The transfer function need 
not be the same for each layer, but for ease of discussion, the 
same symbol A is used for neurons in the hidden layer and in 
the output layer. 

During training of an adaptive neural network, a training 
set of input data is provided, and the weights (and perhaps 
biases) are updated based upon some desired output and 
criterion of goodness. For example, suppose for an adaptive 
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neural network there is a set of input variables {x ls x, . . . , x^} 
and a set of output variables {y 1? y 2 , . . . , y^}. It is convenient 

to define an input vector variable x=(x 1; x, . . . , x^) and an 

5 output vector variable y =(y 1; y 2 , . . . , y^). It is also conve- 
nient to define particular realizations of these vectors, where 

we define input vectors x(i)=(x(i) l5 x(i), . . . , x(i) iV ) and 
output vectors y (i)=(y(i) 1? y(i), . . . , y®^), with the index i 

to denoting a particular realization. For example, {x(i), i=l, 
2, . . . , T} may represent T input training data vectors and 

{ y (i), i=l, 2, . . . , T} may represent the resulting T output 
vectors given by an adaptive neural network for some given 
1 5 set of weights (and also perhaps for some given set of biases) . 
For ease of notation, the dependency of the output vectors on 
the set of weights (and perhaps biases) is not shown. 

During training, for each 7(i) there is a corresponding 

20 desired response vector cf(i). For example, suppose the pat- 
tern recognition function performed by an adaptive neural 
network is to map the input into one of two classes. That is, 
there are two patterns to recognize. For a particular example 
in which there are two output neurons so that the dimension of 

25 is N=2, the desired response may be taken as cf(i)=(10) if 

7(i) belongs to one of the two classes, and cf (i)=(0 1) if 7(i) 
belongs to the other one of the two classes. A criterion of 
goodness may be to find the set of weights (and perhaps 
biases) that minimize the sum of errors e(i) over the training 

set { \(i), i=l, 2, . . . , T}, where e(i)=||y (i)- cf(i)||. 

For arbitrarily dimensioned vectors and desired responses, 
the above-described minimization is a well known math- 
ematical problem, and various mathematical techniques for 
3 5 finding the set of weights (and biases) that sati sfy the criterion 
of goodness are well known. For example, the method of 
steepest descents may be used, which may be used in con- 
junction with the error back-propagation neural network 
learning algorithm, well known in the art of adaptive neural 
40 networks. 

For an adaptive neural network with a sigmoid-logarithmic 
transfer function as discussed previously, some embodiments 
may utilize the back-propagation algorithm for training as 
follows. First, train the adaptive neural network using a sig- 
45 moid transfer function until some criterion of goodness is 
satisfied. For example, some set of training input vectors 

{ "x(i), i=l, 2, . . . , T), desired responses (cf(i), i=l, 2, . . . , T}, 
and threshold A is selected, where the initial set of weights are 
chosen randomly. The back-propagation algorithm is run 
50 until 2 z=1 r e(i)<A is satisfied. Second, use the resulting 
weights as an initial set of weights for another training set 
(which may or may not be different from the first training set), 
but now where the sigmoid-logarithmic transfer function is 
used in the back-propagation algorithm. 

55 Once the adaptive neural network has been trained, it may 
then be operated with static weights (pattern recognition 
mode) to perform pattern recognition. Post processing may be 
applied to the output vector from the output neurons. For 
some embodiments, a winner-take-all module may be 
60 applied, whereby the output neuron with the largest output is 
chosen as the winner. Some embodiments may perform a 
multiple winner-take-all module, whereby the next “highest” 
neuron after the winner is selected, and so on for other neu- 
rons. The final outcome (result) of the adaptive neural net- 
65 work may be represented by some L bit number, where 
L=[log 2 (N)], denoting the selected neuron, where the bracket 
denotes the smallest integer larger than or equal to log 2 (N). 
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Some embodiments may perform signal processing algo- 
rithms other than those described previously. For example, 
note that the sum h=2 I - =1 ' M w f x f may be viewed as an inner 

product of a weight vector w=(w x , w 2 , . . . , w^) with the input 

vector x . The weight vector w may also be referred to as a 
codevector. For some embodiments, a processing neuron 
finds the square of the Euclidian distance, denoted as d, 

between an input vector x and a codevector w, that is, d(x , 

w)=|| x-w|| 2 . The function d may be termed a distortion. It is 
passed on as input to the transfer function. That is, the output 
of the neuron employing a distortion measure is A(d). For 
such embodiments in which the distortion is calculated, the 
winning neuron will have a minimum output. Alternatively, 
the output may be taken as 1/A(d), so that the winning neuron 
will have a maximum output. 

The weights in the weight vector w=(w 1; w 2 , . . . , w^) are 
sometimes also referred to as synapse weights. The process- 
ing of x • w or d( x , w)=|| x - w|| 2 may be considered as part of 
a synapse cell, where the neuron cell involves applying the 
transfer function to the result of the synapse cell. However, 
there is no conceptual difference whether or not the synapse 
cell function is considered part of a neuron cell, or is separated 
out from the neuron cell, although there may be implemen- 
tation differences in realizing the processing in hardware. 
Some embodiments may perform processing other than the 

inner product x\ w or the distortion dfx , w)=||7-w|| 2 , so that 
more generally, some embodiments may pass on to a neuron 

some value f(x , w) where f is some function mapping two 
vectors into a number. 

Furthermore, some embodiments may employ a learning 
function other than a conventional back-propagation algo- 
rithm commonly used in adaptive neural networks. For 
example, some embodiments may perform the following pro- 
cessing operations during the learning mode of an adaptive 
neural network. 

Index the weight vectors (codevectors) as w z where the 
index i refers to the neuron (or synapse cell if that terminology 
is being used). Furthermore, it is useful to add another index 

to w . to designate a particular learning iteration, where w .(t) 
refers to the codevector for neuron i at iteration t. 

Associate with each neuron i a winning frequency f z . It is 
also convenient to index f z according to the iteration index, so 
that f z (t) refers to the winning frequency for neuron i at itera- 
tion index t. Following the convention that t=0 for the first 

iteration, initialize w z (0) by choosing them from a set of 
random (or pseudorandom) numbers. Set f z (0)=l for each i. 

Compute the distortion d z (t) (we have also indexed d 
according to the neuron index and the iteration index) where 

d z (t)=d(w z (t), x z (t)) for each neuron (note that we have also 

indexed the input vector x to refer to neuron i and the itera- 
tion index.) Select the neuron with the smallest distortion and 
set its output, denoted as O z (t), as follows (a value of 1 is 
considered high): O z (t)=l if d z (t)<d J (t), for l^i, j^N, and 
O z (t)=l otherwise (there are N neurons). Update the weight 
vectors (codevectors) with the following frequency- sensitive 

training rule and associated winning frequency: w z (t+l)= 

w,(t)+S(t)0,.(t)[x,(t)-w,(t)]; where S(t)=l/f;.(t) if lSf,(t) 
=f TH , and S(t)=0 otherwise; and f z (t+l)=f z (t)+O z (t). S(t) is the 
frequency-sensitive learning rate, and f ri7 is the upper-thresh- 
old frequency. Notice that only the winning codevector is 
updated. The training rule moves the winning codevector 
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toward the training vector by a fractional amount which 
decreases as the winning frequency increases. If f z (t) is larger 
than f TH , then S(t) is set to zero and no further training will be 
performed for the corresponding neuron. 

5 The above operations are performed for the set of training 
vectors. Use of the upper-threshold frequency may avoid 
codevector under-utilization during the training process for 
an inadequately chosen initial codebook of codevectors. The 
selection of the upper-threshold frequency is heuristic and 
10 depends on source data statistics and the training sequence. 
Empirically, an adequate i TH may be chosen to be two to three 
times larger than the average winning frequency. The initial 
codebooks may be created from a pseudorandom number 
generation function. 

15 A feedforward adaptive neural network is amenable to a 
parallel processing architecture because all of the neurons in 
any one layer may process data concurrently. If the function- 
ality of providing the inner product, distortion function, or 
other types of functions involving the weight vector and input 
20 vector to a neuron is to be separated out from the neuron and 
realized by separate circuits, e.g., the synapses as discussed 
previously, then the previous sentence should be modified to 
indicate that the synapses for a layer also may process data 
concurrently. Furthermore, for multiple layers in a feedfor- 
25 ward adaptive neural network, there may be further concur- 
rency in the sense that one or more layers (of neurons and 
synapses) may be processing in a pipelined fashion. As a 
result, an adaptive neural network is suitable for a VLSI (Very 
Large Scale Integration) circuit that takes advantage of con- 
30 current (parallel) processing. 

The functional blocks of VLSI circuit 401 according to an 
embodiment are illustrated in FIG. 4. The analog voltage 
signals from the differential logarithmic amplifiers (e.g., the 
circuit illustrated in FIG. 1) are provided at input port 402. 
35 These analog signals are provided to sample and hold 404. 
Host Processor 406 provides several functions to VLSI circuit 
401. For example, host processor provides weight vectors to 
sample and hold 408 . These weight vectors may be the weight 
vectors obtained after the adaptive neural network has been 
40 trained, in which case they are provided to synapse weight 
matrix unit 410; or they may be the weight vectors that are 
used for training during a learning mode, in which case they 
are provided to synapse weight matrix 412. During learning, 
a training set is provided to input port 402. The term matrix is 
45 used in the description for functional units 410 and 412 
because the weight vectors may be considered rows (or col- 
umns) of a matrix. 

Note that for the embodiment of FIG. 4, VLSI circuit 
performs much of the learning algorithm, so that the parallel 
50 processing available from functional units 410, 412, and 420 
may be utilized, in which case host process 406 performs 
some non-parallel learning functions and data-flow control. 
For some embodiments, all or some of the functions provided 
by host processor 406 may be integrated on VLSI circuit 401 . 
55 Control lines 413 allow host processor 406 to select the 
sources of input vectors to VLSI circuit 401, and whether 
weight matrix 410 (pattern recognition mode) or 412 (adap- 
tive or learning mode) are used to store the weight vectors . For 
example, for some embodiments, if for control lines 413 “IV” 
60 (mnemonic for “input vector”) is set to a logic 0 (a LOW 
digital signal value), then the input vector is from the host 
processor; whereas if “IV” is set to a logic 1 (a HIGH digital 
signal value), then the input vector is from sensor array 110 , 
that is, the analog output of the differential logarithmic ampli- 
65 fier in circuit 114. If for control lines 413 “WV” (mnemonic 
for “weight vector”) is set to a logic 1 , then weight matrix 410 
is selected to store weight vectors loaded from host processor 
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406 . These weight vectors are those that are obtained after the 
learning algorithm has been performed, so that the adaptive 
neural network is operating in its pattern recognition mode; or 
they may be the desired response vectors d(i) used in a super- 
visory learning algorithm when the neural network is operat- 
ing in its adaptive or (supervisory) learning mode. 

The particular weight vector in functional units 410 and 
412 for a neuron is addressed by address decoders 414 and 
416 , where the particular address is provided by host proces- 
sor 406 by way of vector address bus (lines) 417 . Input neu- 
rons 418 provides the analog input vector to either synapse 
weight matrix 410 or 412 , depending upon whether the adap- 
tive neural network is in a learning mode or a pattern recog- 
nition mode. The latter may be termed an encoding mode, in 
the sense that an input vector is encoded into a recognizable 
class. Functional units 410 or 412 perform the synapse func- 
tion, where for the previously described embodiments may 
involve forming the inner product of the weight vectors with 
input vectors, or calculating the distortions. These results are 
passed to functional unit 420 . 

Functional unit 420 performs the neuron functions dis- 
cussed previously, that is, functional unit 420 applies the 
transfer function to the synapse result. The particular transfer 
function is selectable. For example, the transfer function may 
be a sigmoid function, or a sigmoid- logarithm as discussed 
previously. The outputs of these neurons are provided to 
functional unit 422 , which performs a winner-take-all func- 
tion, or perhaps selects one among the top several neurons. 
This result may be encoded into a binary number, provided at 
output port 424 . 

Sensor array 110 , circuit 114 , and VLSI circuit 401 may be 
integrated on a single die (system-on- chip) for some embodi- 
ments, whereas for other embodiments these components 
may reside on two or more die, or comprise a multi-chip 
module, for example. 

Although the subject matter has been described in lan- 
guage specific to structural features and methodological acts, 
it is to be understood that the subject matter defined in the 
appended claims is not necessarily limited to the specific 
features or acts described above. Rather, the specific features 
and acts described above are disclosed as example forms of 
implementing the claims. Accordingly, various modifications 
may be made to the described embodiments without depart- 
ing from the scope of the invention as claimed below. 

It is to be understood in these letters patent that the meaning 
of “A is connected to B”, where A or B may be, for example, 
a node or device terminal, is that A and B are connected to 
each other so that the voltage potentials of A and B are 
substantially equal to each other. For example, A and B may 
be connected together by an interconnect (transmission line). 
In integrated circuit technology, the interconnect may be 
exceedingly short, comparable to the device dimension itself. 
For example, the gates of two transistors may be connected 
together by polysilicon, or metal interconnect, where the 
length of the polysilicon, or metal interconnect, is compa- 
rable to the gate lengths. As another example, A and B may be 
connected to each other by a switch, such as a transmission 
gate, so that their respective voltage potentials are substan- 
tially equal to each other when the switch is ON. 

It is also to be understood in these letters patent that the 
meaning of “A is coupled to B” is that either A and B are 
connected to each other as described above, or that, although 
A and B may not be connected to each other as described 
above, there is nevertheless a device or circuit that is con- 
nected to both A and B. This device or circuit may include 
active or passive circuit elements, where the passive circuit 
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elements may be distributed or lumped-parameter in nature. 
For example, A may be connected to a circuit element that in 
turn is connected to B. 

It is also to be understood in these letters patent that a 
5 “current source” may mean either a current source or a current 
sink. Similar remarks apply to similar phrases, such as, “to 
source current”. 

It is also to be understood in these letters patent that various 
circuit components and blocks, such as current mirrors, 
10 amplifiers, etc., may include switches so as to be switched in 
or out of a larger circuit, and yet such circuit components and 
blocks may still be considered connected to the larger circuit. 

Throughout the description of the embodiments, various 
mathematical relationships are used to describe relationships 
15 among one or more quantities. For example, a mathematical 
relationship or mathematical transformation may express a 
relationship by which a quantity is derived from one or more 
other quantities by way of various mathematical operations, 
such as addition, subtraction, multiplication, division, etc. Or, 
20 a mathematical relationship may indicate that a quantity is 
larger, smaller, or equal to another quantity. These relation- 
ships and transformations are in practice not satisfied exactly, 
and should therefore be interpreted as “designed for” rela- 
tionships and transformations. One of ordinary skill in the art 
25 may design various working embodiments to satisfy various 
mathematical relationships or transformations, but these rela- 
tionships or transformations can only be met within the tol- 
erances of the technology available to the practitioner. 

Accordingly, in the following claims, it is to be understood 
30 that claimed mathematical relationships or transformations 
can in practice only be met within the tolerances or precision 
of the technology available to the practitioner, and that the 
scope of the claimed subject matter includes those embodi- 
ments that substantially satisfy the mathematical relation- 
35 ships or transformations so claimed. 

What is claimed is: 

1. A system comprising: 

4Q an array of pixels, each pixel comprising: 

a first photodetector; 

a first optical filter having a first passband to pass to the first 
photodetector a first filtered light having a first intensity; 

a second photodetector; 

45 a second optical filter having a second passband different 
from the first passband, the second optical filter to pass 
to the second photodetector a second filtered light hav- 
ing a second intensity; and 

a first circuit coupled to the first and second photodetectors 
to provide a voltage indicative of a logarithm of the ratio 
of the first intensity to the second intensity, wherein the 
array of pixels has M pixels, where M is an integer 
greater than one, and the system further comprises a 
second circuit comprising: 

55 a functional unit to store a set of weight vectors w z , i=l, 
2, . . . , N, where N is an integer greater than one, each 

weight vector w z of dimension M, the functional unit to 
calculate quantities h*, i=l, 2, . . . , N where h z =;f(x ,w z ), 

60 where x is an M dimensional vector of the voltages 
provided by the first circuit of each corresponding pixel, 
and f is a function of two M dimensional vectors; and 

a set of N neuron processors, each neuron processor to 
provide a quantity A(h z ) where A is a selectable transfer 
65 function. 

2 . The system as set forth in claim 1 , wherein the selectable 
transfer function may be selected as: 
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-adn(£(<5 - A)) 
orln + A)) 


(l+exp(-A)) 1 


A< 2 

-2 < A < 0 
0< A<2 
2 < A 
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where a, (3, and 6 are constants. 

3. The system as set forth in claim 1, wherein the function 

f is an inner product function where 0(7, w z )=7*w z .. 

4. The system as set forth in claim 1, wherein the function 

f is a distortion function where 0(x ,w z )=|| x -w z || 2 . 

5. The system as set forth in claim 1, further comprising a 
die, wherein the second circuit and the array of pixels are 15 
integrated on the die. 

6. The system as set forth in claim 1, further comprising: an 
assay array comprising wells; and a lens system to provide an 
optical path from the assay array to the array of pixels. 

7. A system comprising a set of M pixels pixel(i), i=l, 20 
2, . . . , M, where M is an integer greater than one, for each i=l , 

2, . . . , M, pixel(i) comprising: 
a first photodetector(i); 

a first optical filter(i) having a first passband to pass 
through imaged light to the first photodetector(i); 25 

a second photodetector(i); 

a second optical filter(i) having a second passband different 
from the first passband to pass through imaged light to 
the second photodetector(i); 

a first transistor(i) having a drain connected to the first 30 
photodetector(i) and having a gate connected to the drain 
of the first transistor(i); 

a second transistor(i) having a drain connected to the sec- 
ond photodetector(i) and having a gate connected to the 
drain of the second transistor(i); 35 

a third transistor(i) having a gate connected to the gate of 
the first transistor(i) and having a source; and 
a fourth transistor(i) having a gate connected to the gate of 
the second transistor(i), having a source connected to the 
source of the third transistor(i), and having a drain to 40 
provide a voltage x z . 

8. The system as set forth in claim 7, further comprising: 

a functional unit to store a set of weight vectors w z , i=l , 

2, . . . , N, where N is an integer greater than one, each ^ 

weight vector w z of dimension M, the functional unit to 
calculate quantities h z , i=l, 2, . . . , N where h z =;f (7, w z ), 

where x is an M dimensional vector with component i 
equal to the voltage x z , and f is a function of two M 
dimensional vectors; and 

a set of N neuron processors, each neuron processor to 
provide a quantity A(h z ) where A is a selectable transfer 
function. 

9. The system as set forth in claim 7, further comprising a 55 
die, wherein the functional unit, the set of N neuron proces- 
sors, and the set of M pixels arc integrated on the die. 

10. The system as set forth in claim 7, further comprising: 
an assay array comprising wells; and 

a lens system to provide an optical path from the assay 
array to the set of M pixels. 
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1 1 . The system as set forth in claim 7, the third transistor (i) 
having a drain, the system further comprising for each i=l, 
2, . . . , M, pixel(i): 

a first impedance device(i) connected to the drain of the 
third transistor(i); and 

a second impedance device(i) connected to the drain of the 
fourth transistor(i). 

12. The system as set forth in claim 8, wherein the select- 
able transfer function may be selected as: 


m = 


(1 + exp(-A)) _1 
-adn(/?(<5 - A)) 
adn(/?(<5 + A)) 


(1 + exp(-A)) 1 


A< 2 

-2 < A < 0 
0 < A < 2 
2 < A 


where a, p, and 8 are constants. 

13. The system as set forth in claim 8, wherein the function 

f is an inner product function where f( x ,w z )=x *w z . 

14. The system as set forth in claim 8, wherein the function 

f is a distortion function where ;f(7,w z )=||7-w z || 2 . 

15. The system as set forth in claim 11, wherein for each 
i=l, 2, . . . , M: 

the first impedance device(i) is a first resistor(i); and 

the second impedance device(i) is a second resistor(i). 

16. The system as set forth in claim 11, the system further 
comprising for each i=l, 2, . . . , M, pixel(i): 

a tail transistor(i) connected to the sources of the third 
transistor(i) and the fourth transistor(i). 

17. An adaptive neural network, comprising: 

a functional unit to store a set of weight vectors w z , i=l, 
2, . . . , N, where N is an integer greater than one, each 

weight vector w z - of dimension M, where M is an integer 
greater than one, the functional unit to calculate quanti- 
ties h z , i=l, 2, . . . , N where h^j^x ,w z ), where x is an 
M dimensional vector, and f is a function of two M 
dimensional vectors; and 

a set of N neuron processors, each neuron processor to 
provide a quantity A(h z ) where A is a selectable transfer 
function, wherein the selectable transfer function may 
be selected as: 


(l+exp(-A)) 1 
-adn(/?(<S - A)) 
a\n(j3(6 + A)) 

(1 + exp(-A)) -1 


A < 2 

-2 < A < 0 
0 < A < 2 
2 < A 


where a, p, and 8 are constants. 

18. The system as set forth in claim 17, wherein the func- 
tion f is an inner product function where f( x , w z )= x w z . 

19. The system as set forth in claim 17, wherein the func- 
tion f is a distortion function where f( x , w .)=|| x -w .|| 2 . 





